The Casimir effect within scattering theory 
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We review the theory of the Casimir effect using scattering techniques. After years of theoretical 
efforts, this formalism is now largely mastered so that the accuracy of theory-experiment comparisons 
is determined by the level of precision and pertinence of the description of experimental conditions. 
Due to an imperfect knowledge of the optical properties of real mirrors used in the experiment, 
the effect of imperfect reflection remains a source of uncertainty in theory-experiment comparisons. 
For the same reason, the temperature dependence of the Casimir force between dissipative mirrors 
remains a matter of debate. We also emphasize that real mirrors do not obey exactly the assumption 
of specular reflection, which is used in nearly all calculations of material and temperature corrections. 
This difficulty may be solved by using a more general scattering formalism accounting for non- 
specular reflection with wavevectors and field polarizations mixed. This general formalism has 
already been fruitfully used for evaluating the effect of roughness on the Casimir force as well as the 
lateral Casimir force appearing between corrugated surfaces. The commonly used 'proximity force 
approximation' turns out to lead to inaccuracies in the description of these two effects. 



I. INTRODUCTION 



After its prediction in 1948 [1], the Casimir force 
has been observed in a number of 'historic' experiments 
which confirmed its existence and main properties [2-5]. 
With present day technology, a new generation of Casimir 
force measurements has started since nearly a decade ago 
[6-12]. These experiments have reached a good enough 
accuracy to allow for a comparison between theoretical 
predictions and experimental observations which is of 
great interest for various reasons [13-15]. 

The Casimir force is the most accessible effect of vac- 
uum fluctuations in the macroscopic world. As the exis- 
tence of vacuum energy raises difficulties at the interface 
between the theories of quantum and gravitational phe- 
nomena, it is worth testing this effect with the greatest 
care and highest accuracy [16, 17]. A precise knowledge 
of the Casimir force is also a key point in many accurate 
force measurements for distances ranging from nanome- 
ter to millimeter. These experiments are motivated either 
by tests of Newtonian gravity at millimetric distances 
[18-21] or by searches for new weak forces predicted in 
theoretical unification models with nanometric to mil- 
limetric ranges [22-27]. Basically, they aim at putting 
limits on deviations of experimental results from present 
standard theory. As the Casimir force is the dominant 
force between two neutral non-magnetic objects in the 
range of interest, any new force would appear as a differ- 
ence between experimental measurements and theoreti- 
cal expectations of the Casimir force. On a technological 
side, the Casimir force has been shown to become impor- 
tant in the architecture of micro- and nano-oscillators 
(MEMS, NEMS) [28, 29]. In this context, it is extremely 
important to account for the conditions of real experi- 
ments. 

The comparison between theory and experiment 
should take into account the important differences be- 
tween the real experimental conditions and the ideal sit- 



FIG. 1: Original Casimir configuration of two plane parallel 
mirrors a distance L apart. 



nation considered by Casimir. Casimir calculated the 
force between a pair of perfectly smooth, flat and paral- 
lel plates in the limit of zero temperature and perfect re- 
flection (see Fig.l). He found an expression for the force 
^Cas and the corresponding energy -Ecas which only de- 
pend on the distance L, the area A and two fundamental 
constants, the speed of light c and Planck constant H 
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Each transverse dimension of the plates has been sup- 
posed to be much larger than L. Conventions of sign have 
been chosen so that i^cas is positive while -Bcas is neg- 
ative. They correspond to an attractive force (^ O.l^N 
for A = Icm^ and L = 1/im) and a binding energy. 

The fact that the Casimir force (1) only depends on 
fundamental constants and geometrical features is re- 
markable. In particular it is independent of the fine 
structure constant which appears in the expression of the 
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atomic Van der Waals forces. This universality property 
is related to the assumption of perfect reflection used by 
Casimir in his derivation. Perfect mirrors correspond to a 
saturated response to the fields since they reflect 100% of 
the incoming light. This explains why the Casimir effect, 
though it has its microscopic origin in the interaction of 
electrons with electromagnetic fields, does not depend on 
the fine structure constant. 

However, no real mirror can be considered as a perfect 
reflector at all field frequencies. In particular, the most 
precise experiments are performed with metallic mirrors 
which show perfect reflection only at frequencies smaller 
than a characteristic plasma frequency up which depends 
on the properties of conduction electrons in the metal. 
Hence the Casimir force between metal plates can flt the 
ideal Casimir formula (1) only at distances L much larger 
than the plasma wavelength 



For metals used in the recent experiments, this wave- 
length lies in the O.l/um range (107nm for Al and 137nm 
for Cu and Au). At distances smaller than or of the or- 
der of the plasma wavelength, the finite conductivity of 
the metal has a significant effect on the force. The idea 
has been known since a long time [30 32] but a precise 
quantitative investigation of the effect of imperfect re- 
flection has been systematically developed only recently 
[33-36]. As the effect of imperfect reflection is large in 
the most accurate experiments, a precise knowledge of 
its frequency dependence is essential for obtaining an ac- 
curate theoretical prediction of the Casimir force. 

This is also true for other corrections to the ideal 
Casimir formula associated with the experimental con- 
figuration. For experiments at room temperature, the ef- 
fect of thermal field fluctuations, superimposed to that of 
vacuum, affects the Casimir force at distances larger than 
a few microns. Again the idea has been known for a long 
time [37, 38] but a quantitative evaluation taking into 
accoimt the correlation of this effect with that of imper- 
fect reflection has been mastered only recently [39, 40]. 
A number of publications have given rise to contradic- 
tory estimations of the Casimir force between dissipative 
mirrors at non zero temperature [41-45]. Many attempts 
have been made to elucidate the problem by taking into 
account the low-frequency character of the force between 
metallic films [46], the spatial dispersion on electromag- 
netic surface modes [47] or the transverse momentum de- 
pendance of surface impedances [48-50] . Experimentally 
the effect of temperature of the Casimir force has not yet 
be conclusively measured [51]. For a recent analysis of 
this issue see reference [52]. 

Most experiments are performed between a plane and 
a sphere with the force estimation involving a geometry 
correction. Usually the Casimir force in the plane-sphere 
(PS) geometry is calculated using the Proximity Force 
Approximation (PFA). This approximation amounts to 
the addition of force or energy contributions correspond- 



ing to different local inter-plate distances, assuming these 
contributions to be independent. But the Casimir force 
and energy are not additive, so that the PFA cannot be 
exact, although it is often improperly called a theorem. 

In the present review, we consider both the original 
Casimir geometry with perfectly plane and parallel mir- 
rors and the plane-sphere geometry when comparing to 
experiments. The PFA is expected to be valid in the 
plane-sphere geometry, when the sphere radius R is much 
larger than the separation L [53-55] , which is the case for 
all present day experiments, and it will thus be used to 
connect the two geometries. In this case, the force Fps 
between a sphere of radius R and a plane at a distance 
of closest approach L is given in terms of the energy Epp 
for the plane-plane cavity as follows 

Fps = 27ri?:^ , L^R (3) 

Interesting attempts to go beyond this approximation 
concerning the plane-sphere geometry have been made 
recently [56, 57], in the more general context of the con- 
nection between geometry and the Casimir effect [58-60] . 

Another important correction to the Casimir force is 
coming from surface roughness, which is intrinsic to any 
real mirror, with amplitude and spectrum varying de- 
pending on the surface preparation techniques. The de- 
parture from flatness of the metallic plates may also be 
designed, in particular under the form of sinusoidal cor- 
rugation of the plates which produce a measurable lat- 
eral component of the Casimir force [61]. For a long 
time, these roughness or corrugation corrections to the 
Casimir force have been calculated with methods valid 
only in limiting cases [62-67] or by using PFA [68-70]. 
Once again, it is only recently that emphasis has been put 
on the necessity of a more general method for evaluating 
the effect of roughness outside the region of validity of 
PFA with imperfect mirrors at arbitrary distances from 
each other [71]. While the condition L ^ i? is suffi- 
cient for applying the PFA in the plane-sphere geometry, 
more stringent conditions are needed for PFA to hold for 
rough or corrugated surfaces. The surfaces should indeed 
be nearly plane when looked at on a scale comparable 
with the separation L, and this condition is not always 
satisfied in experiments. When PFA is no longer valid, 
the effect of roughness or corrugation can be evaluated 
by using the scattering theory extended to the case of 
non-specular reflection. 

We review the Casimir effect within scattering theory 
and the theory of quantum optical networks. The main 
idea of this derivation is that the Casimir force has its 
origin in a difference of the radiation pressure of vacuum 
flelds between the two mirrors and in the outer free fleld 
vacuum. This vacuum radiation pressure can be written 
as an integral over all modes, each mode being associated 
with reflection amplitudes on the two mirrors. We first 
present formulas written for specular reflection which are 
valid for lossless [72] as well as lossy mirrors [73]. We dis- 
cuss the influence of the mirrors reflection coefficients at 
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FIG. 2: Schematic view of a Fabry-Perot cavity of length L. 



zero and non-zero temperature. We then extend the ap- 
proach to the case of non-specular reflection which mixes 
the field polarizations and transverse wave- vectors. Fi- 
nally we apply the latter approach to the calculation of 
the roughness correction to the Casimir force between 
metallic mirrors [74, 75] and of the lateral component of 
the Casimir force between corrugated plates [76]. 



II. SPECULAR SCATTERING 

Let us first consider the original Casimir geometry with 
perfectly plane and parallel mirrors aligned along the di- 
rections X and y. The two mirrors thus form a Fabry- 
Perot cavity of length L as shown in Fig. 2. We analyze 
the cavity as a composed optical network, and calcu- 
late the fluctuations of the intracavity fields propagat- 
ing along the positive and negative z-axis, Ec and Ecj 
in terms of the fluctuations of the incoming free-space 
fields Ei^" and Ej^^ (the outgoing fields E°"* and E™* are 
also shown). 

The field modes are conveniently characterized by their 
frequency cj, transverse wavevector k with components 
kx , ky in the plane of the mirrors and polarization p. As 
the configuration of Fig. 2 obeys a symmetry with respect 
to time translation as well as transverse space transla- 
tions (along directions x and y), the frequency uj, trans- 
verse vector k = {kx,ky) and polarization p — TE,TM 
are preserved throughout the whole scattering processes 
on a mirror or a cavity. The scattering couples only the 
free vacuum modes which have the same values for the 
preserved quantum numbers and differ by the sign of the 
longitudinal component k^ of the wavevector. We denote 
by (r^ reflection amplitude of the mirror j = 1, 2 

as seen from the inner side of the cavity. This scattering 
amplitude obeys general properties of causality, unitarity 
and high frequency transparency. The additional fluctua- 
tions accompanying losses inside the mirrors are deduced 
from the optical theorem applied to the scattering pro- 



cess which couples the modes of interest and the noise 
modes [77, 78]. 

The loop functions which characterize the optical re- 
sponse of the cavity to an input field play an important 
role in the following 



2ik^L 



(4) 



and are respectively the open-loop and closed-loop 
functions corresponding to one round trip in the cavity. 
The system formed by the mirrors and fields is stable so 
that is an analytic function of frequency uj. Analyt- 
icity is defined with the following physical conditions in 
the complex plane 



oj = iS^ 



SR^ > 



(5) 



kz = in [lo\ 



K[tj]EEWk2~— , SRk[w]>0 



The quantum numbers p and k remain spectator 
throughout the discussion of analyticity. The sum on 
transverse wavevectors may be represented as a sum over 
the eigenvectors kx = 2Trqx/Lx,ky = 2'Kqy/Ly associated 
with virtual quantization boxes along x, y or, at the con- 
tinuum limit LxT Ly oo with A = LxLy, as an integral 
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We then introduce the Airy function defined in clas- 
sical optics as the ratio of energy inside the cavity to 
energy outside the cavity for a given mode 

5n-]-i + {/n-]+--} = ?-4r!P 

/k' 5k depend only on the reflection amplitudes of mir- 
rors as they are seen from the inner side. With these 
definitions we write the Casimir force 
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F = -^EE I ^{-H/kH^] 



+ C.C.} (8) 



p k 

or, equivalently, the Casimir energy 
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(9) 



Equations (8,9) contain the contribution of ordinary 
modes freely propagating outside and inside the cavity 
with Lu > cjkj and kz real. This contribution thus merely 
reflects the intuitive picture of a radiation pressure of 
fluctuations on the mirrors of the cavity [72] with the 
factor 5^ — 1 representing a difference between inner and 
outer sides. Equations (8,9) also include the contribu- 
tion of evanescent waves with uj < c|k| and kz imaginary. 
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Those waves propagate inside the mirrors with an inci- 
dence angle larger than the limit angle and they also exert 
a radiation pressure on the mirrors, due to the frustrated 
reflection phenomenon [73]. Their properties are conve- 
niently described through an analytical continuation of 
those of ordinary waves, using the well defined analytic 
behavior of k and 

Using analyticity properties, we now transform (8) into 
an integral over imaginary frequencies by applying the 
Cauchy theorem on the contour enclosing the quadrant 
^u) > 0, Sjcj > 0. We use high frequency transparency to 
neglect the contribution of large frequencies. This leads 
to the following expression for the Casimir force 



+ C.C.} (10) 



which is now written as an integral over complex fre- 
quencies Lo = i^. In the same way we obtain the Casimir 
energy as a function of imaginary frequencies 



E : 



p 



oc 

J^Jd^ln[l-fim (11) 



Causality and passivity conditions assure that the inte- 
grand In [1 — /o^^ [i^]] is analytical in the upper half space 
of the complex plane > 0. It is thus clear that both 
expressions for force and energy are equivalent. 



A. Finite conductivity correction 

Let us now review the correction to the Casimir force 
coming from the finite conductivity of any material. This 
correction is given by relations (8) or, equivalently, (10), 
as soon as the reflection amplitudes are known. These 
amplitudes are commonly deduced from models of mir- 
rors, in particular bulk mirrors, slabs or layered mirrors, 
the optical response of metallic matter being described 
by some permittivity function. This function may be ei- 
ther a simple description of conduction electrons in terms 
of a plasma or Drude model or a more elaborate repre- 
sentation based upon tabulated optical data. At the end 
of the section, we will discuss the uncertainty in the the- 
oretical evaluation of the Casimir force coming from the 
lack of knowledge of the specific material properties of a 
given mirror as was illustrated in [34] . 

Assuming that the metal plates have a large optical 
thickness, the reflection coefficients correspond to the 
ones of a simple vacuum-bulk interface [80] 



(e (^0-1) + c2k2-ck 
Ve (£(*6-1) + 



c^K^ — cue (i^) 



{e{ii)-l)+c^K-' + CKe{iC) ^"^^^ 

stands for r'^ {i^, m) and e {i^) is the dielectric func- 
tion of the metal evaluated for imaginary frequencies; the 
index k has been dropped. 



Taken together, the relations (10,12) reproduce the Lif- 
shitz expression for the Casimir force [30] . Note that the 
expression was not written in this manner by Lifshitz. To 
our present knowledge, Kats [81] was the first to stress 
that Lifshitz expression could be written in terms of the 
reflection amplitudes (12). We then have to emphasize 
that (10) is much more general than Lifshitz expression 
since it still holds with mirrors characterized by reflec- 
tion amplitudes differing from (12). As an illustration, 
we may consider metallic slabs having a flnite thickness. 

For a given polarization, we denote by Vsi the reflection 
coefficient (12) corresponding to a single vacimm/metal 
interface and we write the reflection amplitude r for the 
slab of finite thickness through a Fabry-Perot formula 



1 - e 
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6 = 



(13) 



This expression has been written directly for imaginary 
frequencies. The parameter 5 represents the optical 
length in the metallic slab and D the physical thickness. 
The single interface expression (12) is recovered in the 
limit of a large optical thickness S ^ 1. With the plasma 
model, this condition just means that the thickness D is 
larger than the plasma wavelength Ap . 

In order to discuss experiments, it may also be worth 
to write the reflection coefficients for multilayer mirrors. 
For example one may consider two-layer mirrors with a 
layer of thickness D oi a, metal A deposited on a large 
slab of metal B in the limit of large thickness as shown 
in Figure 3. The reflection formulas are then obtained as 

A B 



AB 

FIG. 3: Composition of networks : two networks labeled A 
and B are piled up to build up a network AB. 

in [82] but accounting for oblique incidence 



Tab = ta 



1 - rATB 



tAB 



1 - tatb 



(14) 



It reproduces the known results for the simple multilayer 
systems which have already been studied [13]. The com- 
bination of (13) and (14) allows to calculate most of the 
experimental situations precisely. 

In order to assess quantitatively the effect of finite con- 
ductivity, we may in a first approach use the plasma 
model for the metallic dielectric function, with oj-p the 
plasma frequency. 



e{oo) = l- 
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FIG. 4: Reduction of the Casimir force compared to the force 
between perfect mirrors, when the finite conductivity is de- 
scribed by a plasma model (solid line) or a Drude model 
(dashed line) with a ratio ^ equal to 4 x 10~^. The differ- 
ence due to the relaxation parameter has only a small effect 
on the calculation of the Casimir force. The dotted-dashed 
line corresponds to the short distance asymptotic behavior 
(17). 



It is convenient to present the change in the Casimir force 
in terms of a factor "qp which measures the reduction of 
the force with respect to the case of perfect mirrors 



F = r]FFc 



(16) 



Using expressions (12,15) it is possible to obtain the re- 
duction factor defined for the Casimir force through nu- 
merical integrations. 

The result is plotted as the solid line on figure 4, as a 
function of the dimensionless parameter that is the 
ratio between the distance L and the plasma wavelength 
Ap. As expected the Casimir formula is reproduced at 
large distances {r]F — » 1 when L Ap). At distances 
smaller than Ap in contrast, a significant reduction is 
obtained with the asymptotic law of variation read as 
[84, 85] 



i < Af 



r]F = a— a ~ 1.193 (17) 
Ap 



This can be understood as the result of the Coulomb in- 
teraction of surface plasmons at the two vacuum/metal 
interfaces [83, 84]. The generalization of this idea at ar- 
bitrary distances is more subtle since it involves a full 
electromagnetic treatment of the plasmon as well as or- 
dinary photon modes [86] . 

The plasma model cannot provide a fully satisfactory 
description of the optical response of metals, in particu- 
lar because it docs not account for any dissipativc mecha- 
nism. A more realistic representation is the Drude model 
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e(e+7) 



(18) 



This model describes not only the plasma response of 
conduction electrons with wp still interpreted as the 
plasma frequency but also their relaxation, 7 being the 
inverse of the electronic relaxation time. 

The relaxation parameter 7 is much smaller than the 
plasma frequency. For Al, Au, Cu in particidar, the ratio 
j/ujp is of the order of 4 x 10^'\ Hence relaxation affects 
the dielectric constant in a significant manner only at 
frequencies where the latter is much larger than unity. 
In this region, the metallic mirrors behave as a nearly 
perfect reflectors so that, finally, the relaxation docs not 
have a large influence on the Casimir effect at zero tem- 
perature. This qualitative discussion is confirmed by the 
result of numerical integration reported as the dashed 
line on figure 4. With the typical value already given 
for -y/cop, the variation of r]F remains everywhere smaller 
than 2%. 

For metals like Al, Au, Cu, the dielectric constant de- 
parts from the Drude model when interband transitions 
are reached, that is when the photon energy reaches a 
few eV. Hence, a more precise description of the dielec- 
tric constant should be used for evaluating the Casimir 
force in the sub-/xm range. This description relies on one 
hand on the causality relations obeyed by the dielectric 
response function and on another hand on known opti- 
cal data. The reader is referred to [34] for a detailed 
analysis, but we recall here the main argument and some 
important details. Let us first recall that frequencies are 
measured either in eV or in rad/s, using the equivalence 
1 eV = 1.519 X 10^"'' rad/s. An erroneous conversion fac- 
tor 1 eV = 1.537 X 10^'^ rad/s was used in [34], which led 
to a difference in e{i£,) of less than 1% over the relevant 
distance range. In the end of the calculation, this was 
corresponding to a negligible error in the Casimir force 
and energy [88]. 

The values of the complex index of refraction for dif- 
ferent metals, measured through different optical tech- 
niques, are tabulated as a function of frequency in sev- 
eral handbooks [89-91]. Optical data may vary from one 
reference to another, not only because of experimental 
uncertainties but also because of the dispersion of ma- 
terial properties of the analyzed samples. Moreover, the 
available data do not cover a broad enough frequency 
range so that they have to be extrapolated. These prob- 
lems may cause variations of the results obtained for the 
dielectric function e(i^) and, therefore, for the Casimir 
force. 

Figure 5 shows two different plots of e {i^) for Cu as a 
function of imaginary frequency ^. The solid line corre- 
sponds to the first data set with data points taken from 
[89, 90] and extrapolation at low frequency with a Drude 
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FIG. 5: Dielectric function of Cu versus imaginary frequency. 
The solid line corresponds to the first set (optical data given 
by [89, 90], low frequencies extrapolated by a Drude model 
with u>p = 8.97 eV and 7 = 29.5 meV), the dashed line to 
the second set (optical data given by [91], low frequencies 
extrapolated by a Drude model with wp = 7.5 eV and 7 = 
130 me V). 



FIG. 6: Reduction factor rjp for the Casimir force between two 
Cu plates as a function of the plate separation L. The solid 
line corresponds to the first set (optical data given by [89] , low 
frequencies extrapolated by a Drude model with top = 8.97 eV 
and 7 — 29.5 mcV), the dashed line to the second set (optical 
data given by [91], low frequencies extrapolated by a Drude 
model with uip = 7.5 eV and 7 = 130 meV). 



model with parameters o^p = 8.97 eV and 7 = 29.5 meV 
in reasonable agreement with existing knowledge from 
solid state physics. However, as explained in [34], the 
optical data available for Cu do not permit an unam- 
biguous estimation of the two parameters wp and 7 sep- 
arately. Other couples of values can be chosen which are 
also consistent with optical data. To make this point ex- 
plicit, we have drawn a second plot on figure 5 (dashed 
line) with data taken from [91] and the low frequency in- 
terpolation given by a Drude model with ujp = 7.5 eV and 
7 = 130 meV. These values lead to a dielectric function 
e(i^) smaller than in the first data set over the whole 
frequency range, but especially at low frequencies. An 
estimation of the uncertainties associated with this im- 
perfect knowledge of optical data can be drawn from the 
computation of the Casimir force in these two cases. 

Figure 6 shows the reduction factor r/p for the Casimir 
force between two Cu plates as a function of the plate 
separation L for the two sets of optical data. The two cor- 
responding curves have similar dependance on the plate 
separation but the absolute values are shifted from one 
curve to the other. At a separation of 100 nm the dif- 
ference can be as large as 5%. As the plasma frequency 
is basically the frequency above which the mirrors re- 
flectivity diminishes considerably, the Drude parameters 
of the first set {ujp = 8.97 eV and 7 — 29.5 meV) give 
a larger Casimir force than the second set, where the 
plasma frequency is lower {wp = 7.5 eV and 7 = 130 
mcV). A detailed analysis of this uncertainty has been 
recently reported [92]. 

Let us emphasize that the problem here is neither due 
to a lack of precision of the calculations nor to inaccu- 
racies in experiments. The problem is that calculations 



and experiments may consider physical systems with dif- 
ferent optical properties. Material properties of mirrors 
indeed vary considerably as a function of external pa- 
rameters and preparation procedure [92]. This difficulty 
could be solved by measuring the reflection amplitudes 
of the mirrors used in the experiment and then inserting 
these informations in the formula giving the predicted 
Casimir force. In order to suppress the uncertainty as- 
sociated with the extrapolation procedure, it would be 
necessary to measure the reflection amplitudes down to 
frequencies of the order of 1 meV, if the aim is to calcu- 
late the Casimir force in the distance range from 100 nm 
to a few /um. 



B. Temperature correction 

The Casimir force between dissipative metallic mirrors 
at non zero temperature has given rise to contradictory 
claims which have raised doubts about the theoretical ex- 
pression of the force. In order to contribute to the resolu- 
tion of this difficulty, we now review briefly the derivation 
of the force from basic principles of the quantum theory 
of lossy optical cavities at non zero temperature. We 
obtain an expression which is valid for arbitrary mirrors, 
including dissipative ones, characterized by frequency de- 
pendent reflection amplitudes. This expressions coincides 
with the usual Lifshitz expression when the plasma model 
is used to describe the mirrors material properties, but it 
differs when the Drude model is applied. The difference 
can be traced back to the validity of Poisson summation 
formula [40]. 

To discuss the effect of finite temperature we use a 
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theorem which gives the comimitators of the intracavity 
fields as the product of those well known for fields out- 
side the cavity by the Airy function. This theorem was 
demonstrated with an increasing range of validity in [72] , 
[79] and [73] . It is true regardless of whether the mirrors 
are lossy or not. Since it does not depend on the state 
of the field, it can be used for thermal as well as vacuum 
fluctuations. Assuming thermal equilibrium, the theo- 
rem leads to the expression of the field anticommutators, 
i.e. the field fiuctuations. Note that thermal equilibrium 
has to be assumed for the whole system, which means 
that input fields as well as fluctuations associated with 
electrons, phonons and any loss mechanism inside the 
mirrors correspond to the same temperature T, whatever 
their microscopic origin may be. If parts of the system 
correspond to different temperatures, completely differ- 
ent results are obtained [94, 95] . 

The anticommutators of intracavity fields are given by 
those known for fields outside the cavity multiplied by 
the Airy function. Hence, the expression written in [73] 
for a null temperature is only modified through the ap- 
pearance of a thermal factor in the integrand 
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may also be read as 



COth ( I , WT = (19) 



Using as before analyticity properties, we transform 
(19) into an integral over imaginary frequencies giving 
the following expression for the Casimir force 
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4i{x) = 2 dCcos(^x)</.£[^] 



lim Kk H + V] /k H + ^] 
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We have introduced the common summation convention 



(22) 



n n=l 

The function 0^ is well defined almost everywhere, the 
only possible exception being the point ^ = where the 
limit r? — > 0+ may be ill defined for mirrors described 
by dissipative optical models [44]. Since this is a do- 
main of null measure, the cosine Fourier transform 0^ of 
is well defined everywhere and the expression (21) of 
the Casimir force is valid for arbitrary mirrors, includ- 
ing dissipative ones. Note that the term n = in (21) 
corresponds exactly to the contribution of vacuum fluctu- 
ations, or to the zero temperature limit, while the terms 
n > 1 give the corrections associated with thermal fields. 

We come back to the derivation of the Lifshitz formula 
[30] , often used as the standard expression of the Casimir 
force. This formula is directly related to the decompo- 
sition of the coth function into elementary fractions cor- 
responding to the Matsubara poles = imiVT- If we 
assume furthermore that the function ^ is a sufficiently 
smooth test function, in the sense defined by the theory 
of distributions, we deduce that the expression (20) can 
also be read 



It is now written as an integral over complex frequencies 
ui = + T] close to the imaginary axis, with the small 
positive real number 77 0"*" maintaining the Matsubara 
poles ojm = imojT of c [u>] outside the contour used to 
apply the Cauchy theorem. Up to this point, the present 
derivation is similar to Lifshitz' demonstration [30] while 
being valid for arbitrary reflection amplitudes. The next 
steps in Lifshitz' derivation, scrutinized in [40], may raise 
difficulties for dissipative mirrors. Let us briefly recall the 
main arguments of [40]. 

We may flrst write a series expansion of the Casimir 
force (19) based upon the expansion of the function 

coth (^^7) ^^^^ ^ series of exponentials cxp ^— ^""^ ^ (see 

also [39]). This expansion obeys the mathematical crite- 
rion of uniform convergence so that, when it is inserted in 

(19), the order of the summation over n and integration 
over ^ may be exchanged. It follows that the force (19) 



^-f = ^EEE'^.M (23) 



This is the generalization of the Lifshitz' formula [30] to 
the case of arbitrary reflection amplitudes. It is a discrete 
sum over Matsubara poles with the primed summation 
symbol having the definition (22). This formula is known 
to lead to the correct result in the case of dielectric mir- 
rors (for which it was derived in [30]), for perfect mirrors 
[37, 38] and also for metallic mirrors described by the 
lossless plasma model [39]. 

However its applicability to arbitrary mirrors remains 
a matter of controversy [44] . The point is that the deriva- 
tion of the Lifshitz' formula (23) requires that the func- 
tion be a sufficiently smooth test function, in the sense 
defined by the theory of distributions. Whether or not 
this is the case at ^ = for calculated from dissi- 
pative optical models constitutes the central question of 
the controversy on the value of the term p=TE, m = 
in Lifshitz' sum [41-45]. Let us repeat that (21) is still 
a mathematically valid expression of the Casimir force 
even when cf)^ is ill defined for ^ in a domain of null mea- 
sure. The question of validity of Lifshitz' formula (23) 
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may also be phrased in terms of applicability of the Pois- 
son summation formula [93]. This applicability depends 
on a smoothness condition which is met for dielectric 
mirrors, for perfect mirrors and for mirrors described by 
the plasma model and this explains why Lifshitz' formula 
(23) may be used as well as (21) in these cases [40]. 

In order to solve this controversy, it is crucial to im- 
prove our knowledge of the reflection amplitudes at low 
frequencies. As already discussed, the best manner to 
do that is to measure these amplitudes on the mirrors 
used in the experiment at frequencies as low as possi- 
ble. Although the theoretical question of a good model- 
ing of mirrors at low frequencies is certainly of interest 
and needs to be answered, the crucial point for a reliable 
theory-experiment comparison is the necessity of assess- 
ing the real behavior of the mirrors used in the experi- 
ments. 



III. NON SPECULAR SCATTERING 

We will now present a more general formalism to cal- 
culate the Casimir force and energy which takes into ac- 
count non-specular reflection by the plates. Non-specular 
reflection is of course the generic reflection process on any 
mirror while specular reflection is an idealization. 

In order to introduce the more general formula, let us 
first rewrite expression (11) of the Casimir energy be- 
tween two flat plates as a sum over modes labeled by the 
$, and m = k, j3 

oo 


A£ [i^] = In (1 - rir2e-2«^) (24) 

This can be interpreted as the energy stored inside the 
cavity during the scattering process. It is expressed in 
terms of the phase shifts [iS] acquired by the field 
modes upon scattering on the cavity. These phase shifts 
arc deduced from the 5-matrix of the cavity [72] in such 
a manner that the Casimir energy is simply equal to the 
logarithm of the determinant of the S— matrix. Using 
the techniques of quantum field theory [59] , this can also 
be written as the trace of matrix, here diagonal, defined 
on these modes 

oo 

Esp = hJ^Yl 1^(1 - ^1^26-^'^^^) \m) (25) 

Here, ri and r2 arc diagonal matrices which contain 
as their diagonal elements the specular reflection ampli- 
tudes, as they are seen from fields inside the cavity, 

(to] r, \m') = Sm,m' n it m) 

=47r25(k-k')5p,p' (26) 



while K is a matrix diagonal over the same modes 

(m| K \m') = ^/K, + e^m,m' (27) 

It is now easy to write down a more general formula 
of the Casimir energy for the case of stationary but non- 
specular scattering 

oo 

^nsp ^hj ^Tr In (1 - Tlie-'^^Tlse-'^^) (28) 



The two matrices TZi and 1^2 are no longer diagonal on 
plane waves since they describe non specular reflection 
on the two mirrors. The propagation factors remain di- 
agonal on plane waves. Note that the matrices appearing 
in (28) no longer commute with each other. As a con- 
sequence, the two propagation matrices in (28) can be 
moved through circular permutations in the product but 
not adjoined to each other. 

Formula (28) has already been used to evaluate the 
effect of roughness [74, 75] or corrugation [76] of the mir- 
rors on the Casimir force. To this aim, it was expanded 
at second order in the profiles of the mirrors, with the op- 
tical response of the bulk metals described by the plasma 
model. The non specular reflection amplitudes were then 
deduced from techniques developed for treating rough 
plates [96, 97]. The condition of validity of this expan- 
sion is that the roughness or corrugation amplitude is the 
smallest of length scales involved in the problem. In this 
regime, it was possible to investigate various domains for 
the roughness or corrugation wavelength and thus to in- 
vestigate the effect of roughness or corrugation outside 
as well as inside the range of validity of the Proximity 
Force Approximation. 

We may again emphasize at this point that the for- 
mula (28) has a wider range of validity than used in 
those applications. It can in principle describe mirrors 
with nanostructured surfaces corresponding to large am- 
plitudes which cannot be treated as a small perturba- 
tion. It can as well deal with more complicated optical 
responses which are described neither by a plasma nor by 
a Drude model. As was extensively discussed above for 
the case of specular reflection, the formula (28) remains 
valid for arbitrary mirrors, the only problem being to 
obtain the precise form of the reflection matrices to be 
inserted into it. 



A. Influence of surface roughness 

Let us now recall how the non-specular scattering for- 
mula (28) can be used to calculate the effect of rough- 
ness on the Casimir force. Taking this effect into account 
simultaneously with that of finite conductivity is essen- 
tial, because both of them are important at short dis- 
tances. In order to analyze the roughness effect between 
two metallic plates, we will describe the optical proper- 
ties of the mirrors by the plasma model. The values for 
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z = ^ 
FIG. 7: Magnified detail of the internal surface of mirror Ml. 



the plasma wavelength, the mirror separation and the 
roughness correlation length will be arbitrary with re- 
spect to each other, the roughness amplitude remaining 
the smallest length scale for perturbation theory to hold. 
We will review some simple analytical expressions for sev- 
eral limiting cases, as well as numerical results allowing 
one for a reliable calculation of the roughness correction 
in real experiments [74, 75]. 

In a plane-plane geometry, the surface profiles are de- 
fined by the functions hi{x,y) {i = 1,2) giving the local 
heights with respect to the mean separation L along the 
z direction as shown in Fig. 7. These functions are de- 
fined so that they have zero averages. We consider the 
case of stochastic roughness characterized by spectra 



,-(k) 



(ji^v)h,m , z = l,2 (29) 



We suppose the surface A of the plates to contain many 
correlation areas, which allows us to take ensemble or 
surface averages interchangeably. The two plates are con- 
sidered to be made of the same metal and the crossed cor- 
relation between their profiles is neglected ((Ti2(k) = 0). 

We obtain the following variation of the Casimir energy 
Epp up to second order in the perturbations hi [71] 



pp 



d^k 



a(k)a(k) 



47r2 

O-(k) = (7ll (k) + (722 (k) 



(30) 



With our assumptions, the spectrum cr(k) fully charac- 
terizes the roughness of the two plates. The correlation 
length is defined as the inverse of its width. The 
response function Gr(k) then describes the spectral sen- 
sitivity to roughness of the Casimir effect. Symmetry re- 
quires that it only depends on k — |k|. The dependance 
of Gr on k reflects that not only the roughness amplitude 
but also its spectrum plays a role in diffraction on rough 
surfaces [96, 97]. The formula (30) has been obtained for 
the energy in the plane-plane configuration but it also 



determines the force correction JFps in the plane-sphere 
configuration since the PFA is still used for describing 
the weak curvature of the sphere (see below). 

We now focus our attention on the validity of PFA 
for treating the effect of roughness and notice that this 
validity only holds at the limit of smooth surface profiles 
A; — > 0. In fact, the following identity is obeyed by our 
result [74], for arbitrary values of L and Ap, 



Gr (fc ^ 0) = 



(31) 



where the derivative is taken with respect to the plate 
separation L. If we now suppose that the roughness 
spectrum (T(fc) is included inside the PFA sector where 
Gr (fc) — Gr (0), Gr may be replaced by Gr(0) and fac- 
tored out of the integral (30) thus leading to the PFA 
expression [71] 
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pp 



-^p(i) 2 

2 ° 

d^k 



(32) 



- -^a(k) ^ {h{ + hi) 



Air 



In this PFA limit, the correction depends only on the 
variance of the roughness profiles, that is also the in- 
tegral of the roughness spectrum. 

In the general case in contrast, the sensitivity to rough- 
ness depends on the wavevector k. This key point is 
emphasized by introducing a new function Pr(fc) which 
measures the deviation from the PFA [71] 



p,(fc) = 



G,(fc) 
Gr(0)' 



(33) 



This function is plotted on Fig. 8 for several values of L. 
As for all numerical examples considered below, we take 
Ap — 137nm which corresponds to gold covered plates. 

The ratio Pv{k) is almost everywhere larger than unity, 
which means that the PFA systematically underestimates 
the roughness correction. The inlet shows Pr(fc) for small 
values of k where the PFA is a good approximation. To 
give a number illustrating the deviation from the PFA, 
we find pr ~ 1.6 for L = 200nm and k = 0.02nm^^, which 
means that the exact correction is 60% larger than the 
PFA result for this intermediate separation and a typical 
roughness wavelength 27r/fc ~ 300nm. 

Fig. 8 indicates that Pr(fc) grows linearly for large val- 
ues of k. This is a general prediction of our full calcula- 
tions [74] , for arbitrary values of L and Ap , 



Pi-{k) ~ a,- k for 



27r 1 



(34) 



The dimensionless parameter ar/i depends on 2'kL/Xp 
only, and its expression is given by equation (8) in [74]. 
In Fig. 9, we plot the coefficient ar as a function of L 
with Ap = 137nm. In the limit of short distances, we 
recover the expression which was drawn in [71] from older 
calculations [64] (after a correction by a global factor 2) 

Ap 



0.4492L for fc"^ < L < 



27r 



(35) 
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FIG. 8; Variation of versus fc for L=50, 100, 200, 400 nm 
(from bottom to top cmve). 
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FIG. 9: Variation of the coefficient versus L. The ana- 
lytical result for <C L <C Ap is shown as the dotted line 
and for fe"'^ <C Ap <C i as the dashed line. A comparison 
between this second result (dashed straight line) and the ex- 
act pi{k) (solid line) is shown in the inlet for L = 2fj,m. The 
analytical result pr — Lk/3 predicted by the model of perfect 
reflectors (dotted line) is valid only in the intermediate range 
Ap < < L. 



In the opposite limit of large distances, the coefficient ar 
is found to saturate [74] 

a, = 11^ for (36) 

15 ZTT zw 

It is interesting to note that this result differs from the 
long distance behavior which was drawn in [71] from the 



reanalysis of calculations of the effect of sinusoidal cor- 
rugations on perfectly reflecting plates [66]. Perfect re- 
flectors indeed correspond to the limiting case where Ap 
rather than l//c is the shortest length scale. The follow- 
ing result is obtained in this case [74], which effectively 
flts that of [66], 

p,. = ^Lk for Ap<fc"^<L (37) 

The long-distance behavior is thus given by (36) when 
1/fc < Ap <C L but by (37) when Ap < 1/fc < L. 
The cross-over between these two regimes is shown in 
the inlet of Fig. 9, where we plot p^- as a function of k 
for L — 2/im. The failure of the perfect reflection model 
for 1/fc ^ Ap has been given an interpretation in [74]: 
it results from the fact that not only the incoming field 
mode but also the outgoing one have to see the mirror as 
perfectly reflecting for formula (37) to be valid. 

These numerical results can be used to assess the ac- 
curacy of the PFA applied to the problem of roughness. 
PFA is indeed recovered at the limit of very smooth sur- 
face profiles and the deviation from PFA given by our re- 
sults as soon as the roughness wavevector goes out of this 
limit. The mirrors used in a given experiment have a spe- 
cific roughness spectrum which can, and in our opinion 
must be, measured when the experiments are performed. 
The integral (30) then leads to a reliable prediction for 
the roughness correction, as soon as the spectral sensi- 
tivity G'r(fc) and the real spectrum (T(fc) are inserted into 
it. 



B. Lateral Casimir force component 

The spectral sensitivity G'r(fc) involved in the calcula- 
tion of the roughness correction can be considered as a 
further prediction of Quantum ElectroDynamics, besides 
the more commonly studied mean Casimir force, so that 
the comparison of its theoretical expectation with exper- 
iments is an interesting prospect. But this comparison 
can hardly rely on the roughness correction (30) which 
remains in any case a small variation of the longitudinal 
Casimir effect. A more stringent test can be performed 
by studying the lateral component of the Casimir force 
which arises between corrugated surfaces. This lateral 
Casimir force would indeed vanish in the absence of sur- 
face corrugation so that the expression of the spectral 
sensitivity will thus appear directly as a factor in front of 
the lateral Casimir force. For reasons which will become 
clear below, the spectral sensitivity involved in the calcu- 
lation of corrugation effect is a different fimction Gc(fc). 

Nice experiments have shown the lateral Casimir 
force to be measurable at separations of a few hundred 
nanometers [61], that is of the same order of magnitude 
as the plasma wavelength Ap. It follows that these exper- 
iments can neither be analyzed by assuming the mirrors 
to be perfect reflectors [67], nor by using the opposite 
limit of plasmon interaction [64]. It is no more possible 
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FIG. 10: Surface profiles considered for the lateral component 
of the Casimir force. Both surfaces have a sinusoidal corruga- 
tion with ai and 02 being the corrugation amplitudes, b the 
mismatch between the two sinusoidal functions. 



to use the PFA if we want to be able to treat arbitrary 
values of the ratio of the corrugation wavelength Ac to 
the interplate distance L. This is why we emphasize the 
results drawn from the non-specular scattering formula 
(28) which can be used for calculating the lateral Casimir 
force for arbitrary relative values of Ap, Ac and L. The 
only drawback of this calculation is that it is restricted 
to small enough corrugation amplitudes, since the latter 
have to remain the smallest length scale for perturbation 
theory to hold. But the lateral force is known to be ex- 
perimentally accessible in this regime. Again we model 
the optical response of the metallic plates by the plasma 
model. 

The surface profiles of the corrugated plates are de- 
fined by two functions hi{r), with r = (x, y) is the lateral 
position along the surfaces of the plates while i = 1,2 
labels the two plates. As in experiments [61], we con- 
sider the simple case of uniaxial sinusoidal corrugations 
imprinted on the two plates (see Fig. 10) along the same 
direction, say the y direction, and with the same wavevec- 
tor k = 2tt/\c 

hi = fli cos {kx) , h2 — a2 cos {k{x + b)) (38) 

Both profiles hi and /12 have zero spatial averages and 
they are counted as positive when they correspond to 
local length decreases below the mean value L. 

For the purpose of the calculation of the lateral Casimir 
force, the non-specular reflection matrix TZj have to be 
developed up to the first order in the deviations hj from 
flatness of the two plates. They are thus written as 
the sum of a zero-th order contribution identifying to 
the specular reflection amplitude and of a first-order 
contribution proportional to the Fourier component at 
wavevector (k — k') of the surface profiles, this Fourier 
component being able to induce a scattering of the field 
modes from the wavevector k to k' [76]. The correc- 
tion of the Casimir energy SEpp induced by the corru- 
gations arises at second order in the corrugations, with 
crossed terms of the form ai 02 which have the ability 
to induce lateral forces. In other words, the corruga- 
tion sensitivity function Gc{k) obtained below depends 
on the crossed correlation between the profiles of the two 
plates, in contrast to the function Gr(fc) calculated above 
for describing the roughness spectral sensitivity. The lat- 
ter were depending on terms quadratic in hi or /12, and 



their evaluation required that second order non specular 
scattering be properly taken into account. Here, first or- 
der non specular amplitudes evaluated on both plates are 
sufficient. 

The result of the calculation is read as a second-order 
correction induced by corrugations 



SE 



pp 



d^k 

(27r)2A 



Gc(k)i?i(k)i/2(-k) (39) 



with the function Gc(k) given by equation (3) in [76]. 
For isotropic media, symmetry requires Gc (k) to depend 
only on the modulus of the wavevector k = |k|. We may 
also assume for simplicity that the two plates are made 
of the same metallic medium. The energy correction thus 
depends on the lateral mismatch b between the corruga- 
tions of the two plates, which is the cause for the lateral 
force to arise. Replacing the ill-defined (27r)2(5(2) (0) by 
the area A of the plates, we derive from (39) 



SEpp = aia2 cos{kb)Gc{k) 



(40) 



Once again, the result of the PFA is recovered from 
equation (40) as the limiting case /c — > 0, that is also for 
long corrugation wavelengths. This corresponds to nearly 
plane surfaces where the Casimir energy can be obtained 
from the energy .Epp calculated between perfectly plane 
plates by averaging the 'local' distance C = L — hi — h2 
over the surface of the plates. Expanding at second order 
in the corrugation amplitudes and disregarding squared 
terms in and because they cannot produce a lateral 
dependence, we thus recover expression (40) with Gc(fc) 
replaced, for small values of k or equivalently large values 
of Ac, by Gc(0) given by (compare with (31)) 



Ge(fc ^ 0) = 



E^p 



(41) 



This property is ensured, for any specific model of the 
material medium, by the fact that Gc is given by the 
specular limit of non specular refiection amplitudes [76] 
for fc — > 0. For arbitrary values of fc, the deviation from 
PFA is then described by the ratio 



Pc{k) 



Gc(fc) 
Gc(0) 



(42) 



In the following, we discuss explicit expressions of this 
ratio pc given by its general expression (eq. (3) in [76]). 
For the numerical examples, we take Ap = 137nm, corre- 
sponding to gold covered plates. The result is plotted 
on Fig. 11 as a function of k, for different values of the 
distance L. 

For example for a distance L = 50nm, the Proxim- 
ity Force Approximation is correct in the range k < 



O.Olnm^ 



Ac > 628nm) covered by the plot in 



Fig. 11. However, for typical separations of lOOnm or 
larger, pc drops significantly below its PFA value of unity. 
A more detailed discussion can be found in [76] . 

For still larger values of kL, the functions Gc(fc) and 
Pcik) decay exponentially to zero. If we also assume that 
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FIG. 11: Variation of pc versus k with Ap = 137nm and for 
L — 50nm (dotted line), lOOnm (dash-dotted line), 200nm 
(solid line) and 400nm (dashed line). 



C. Comparison to experiments in a plane-sphere 
configuration 

In order to compare the theoretical expression of the 
lateral Casimir force to experiments, we have to con- 
sider the plane-sphere (PS) geometry [61] rather than the 
plane-plane (PP) one. As R :s> L, we use the PFA to con- 
nect the two geometries. Any interplay between curva- 
ture and corrugation is avoided provided that RL ^ Xq . 
These two conditions are met in the experiment reported 
in [61], where R — 100/xm, Ac — 1.2/im and L ~ 200nm. 

We thus obtain the energy correction 6Eps between 
the sphere and a plane at a distance of closest approach 
L as an integral of the energy correction 5Epp in the PP 
geometry 

6EpsiL,b)^ / 6Epp{L',b) (43) 




0.00 0.01 0.02 0.03 0.04 0.05 

FIG. 12: Lateral force amplitude for the plane-sphere setup, 
as a function of k, with figures taken from [61]. The exper- 
imental value k — 0.0052nm~^ is indicated by the vertical 
dashed line. 



kXp 3> 1, we find Gc{k) — ack exp(— fcL) where the 
parameter ac now depends on Ap and L only. This is 
in striking contrast with the behavior of the response 
function for stochastic roughness, which grows linearly 
with k for large k due to the contribution of the second- 
order reflection coefficients [75]. These coefficients do 
not contribute to the second-order lateral effect, which 
is related to two first-order non-specular reflections at 
different plates, separated by a one-way propagation with 
a modified momentum of the order of k. The resulting 
propagation factor is, in the large-fc limit, exp(— kL) w 
exp(— fcL), thus explaining the exponential behavior. 



Then the lateral force is deduced by varying the energy 
correction (43) with respect to the lateral mismatch b be- 
tween the two corrugations. Simple manipulations then 
lead to the lateral Casimir force in the PS geometry 



Filat 
PS 



27raia2 



kRsm{kb) I dL'Gc{k,L') 



(44) 



The force attains a maximal amplitude for sin(A:&) = ±1, 
which is easily evaluated in the PFA regime fc ^ where 
Gc{k) does not depend on fc, so that Fp^' scales as k. As 
k increases, the amplitude increases at a slower rate and 
then starts to decrease due to the exponential decay of 
Gc(fc). For a given value of the separation i, the lateral 
force reaches an optimum for a corrugation wavelength 
such that kL is of order of unity, which generalizes the 
result obtained for perfect reflectors in [66]. In Fig. 12, 
we plot the force Fpg* (for sin(A:6) = 1) as a function of 
k, with figures taken from the experiment of Rcf. [61]. 
We also use the values ai = 59nm and 02 = 8nm of the 
amplitudes for measuring the force as in [61], reminding 
however that our calculations are valid in the perturba- 
tive limit ai, a2 — » 0. 

The plot clearly shows the linear growth for small k 
as well as the exponential decay for large k. The max- 
imum force \s &i k — 0.009nm~^ so that kL ~ 2. The 
experimental value k = 0.0052nm~^ is indicated by the 
dashed line in Fig. 12, and the force obtained as 0.20pN, 
well below the PFA result, indicated by the straight line 
and corresponding to a force of 0.28pN. 

Such a variation in the lateral Casimir force should in 
principle be measurable in an experiment. This could 
lead to the first unambiguous evidence of the limited va- 
lidity of the Proximity Force Approximation, that is also 
to the first observation of a non trivial effect of geometry 
on the Casimir force. 
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IV. CONCLUSION 

In this review we have described the theory of the 
Casimir effect using the techniques of scattering theory. 
We have recalled how this formalism allows ones to take 
into account the real conditions under which Casimir 
force measurements are performed. 

In particular, the finite conductivity effect can be 
treated in a very precise manner, which is a necessity for 
a reliable theory-experiment comparison. There however 
remain inaccuracies in this comparison if the reflection 
amplitudes are drawn from optical models, because of 
the intrinsic dispersion of optical properties of samples 
fabricated by different techniques. We have emphasized 
that these inaccuracies could be circiimvented by mea- 
suring these reflection amplitudes rather than modeling 
them. 

We have then presented the scattering formulation of 
the Casimir force at non zero temperature. This formu- 
lation clears out the doubt on the expression of the force 
while again requiring to have at one's disposal reflection 
amplitudes representing the real properties of the mirrors 
used in the experiments. Let us at this point emphasize 
that the effect of temperature has not been unambigu- 
ously proven in experiments, and that its observation is 
one of the most urgent challenges of experimental re- 



search in the domain. An interesting possibility would 
be to perform accurate measurements of the force at dis- 
tances larger than a few microns, for example by using 
torsional balances [98] . 

In the second part of the paper, we have presented 
a more general scattering formalism which takes into ac- 
count non-specular reflection. We have also discussed the 
application of this formalism for the calculation of the 
roughness correction to the longitudinal Casimir force as 
well as of the lateral component of the Casimir force aris- 
ing between corrugated surfaces. We have argued that 
the spectral sensitivity functions appearing in these ex- 
pressions have to be considered as a new prediction of 
Quantum ElectroDynamics, which differ from the more 
commonly studied mean Casimir force as soon as one 
goes out of the domain of validity of the PFA. This new 
test seems to be experimentally feasible and constitutes 
another challenge of great interest to be faced in the near 
future. 
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